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Abstract 

Liquids with quasi - chemical bonding between molecules are described in terms of vertex 
model. It is shown that this bonding results in liquid - liquid phase transition, which takes 
place between phases with different mean density of intermolecular bonds. The transition 
may be suggested to be a universal phenomena for those liquids. 
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1 Introduction. 

The idea that bonding in molecular liquids may result in existence of two structural states, 
was first suggested for water 0-0 • At low temperatures, supercooled water demonstrates two 
possible local arrangements of molecules at some intermediate length scale, high density (HD) 
and low density (LD) ones. This results in a competition between these two states, so the 
supercooled water may be treated as a binary mixture of HD and LD liquids The liquid - 
liquid phase transition 0l2j in supercooled water is the ultimate manifestation of the HD and LD 
local states. Now, the water HD and LD local structures are well established both in experiment 
and in computer simulations (see,for example, [B]). Intermolecular coupling by hydrogen bonds 
occurs to be a mechanism for HD and LD local states formation [Sj, original idea was published 
in 0. 

Recently, the unexpected physical behavior was found in Benzene, Quinoline, and other 
complex fluids [Zj-^l]. Namely, the anomalies (discontinuities, kinks, peaks) on temperature 
dependencies of a variety of physical properties were observed. These anomalies may be treated 
as a manifestation of changes in the local structure and may be observed, under special experi- 
mental conditions, in many equilibrium and supercooled liquids. The molecules of those liquids 
demonstrate some common features. For considerations related to the structure of the liquid, 
these molecules may be treated as rigid bodies without internal degrees of freedom. The rela- 
tively rigid electronic configurations of the molecules are characterized by a significant separation 
of positive and negative electric charges at the molecular length scale. At very short distance 
between molecules, strong repulsion determines the interaction, and at the only slightly larger 
length-scale, attraction become important. At this distance, the electric charge distribution may 
be characterized in terms of rigid electrical multipoles. Interactions of these multipoles, together 
with the very short-range repulsion, determine the set of mutual positions and orientations of two 
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neighboring molecules corresponding to potential energy minima, which provide intermolecular 
bonds. 

Thus, the intermolecular coupling seems to be a mechanism for two or more structures forma- 
tion at some intermediate length scale. In present paper, we apply known method of description 
of the polymer solution to the statistics of couplings in molecular liquids. The structure of paper 
is as follows. 

• In section 2, we give a brief introduction into formalism applied. 

• Then, in section 3, we analyze the mean - field solution of the model, which demonstrates 
liquid - liquid phase transition. 

• In section 4, exact solution of the model on the Bethe lattice is analyzed. It is shown that 
the phase transition takes place in this case too. 

• In Conclusion, general evaluation of the model is given. 

The aim of the paper is to demonstrate that local structure separation is rather universal phe- 
nomena for liquids with hydrogen - like bonding. 



Statistical description of condensed substances with the possibility of association and aggregation 
of molecules has been intensively developed for a long time. The theory of associated solutions 
[To] , the Flory theory of polymer solutions [HJ], are the milestones on this way. The idea of 
modelling of chemically aggregated systems in terms of scalar order parameter has a long history 
also. Probably I.Lifshitz was first who offered this approach Then, Flory mean - field theory 
[TH] and spin - polymer analogy [IB] were applied to describe polymer systems. In the present 
paper, we consider the model which is a step on this way. 

Quasi - chemical bonding between molecules has been characterized by two main features. 
First, the bond is short - range and directed (i.e., may be drawn as the vector connecting two 
neighboring molecules). Second, the bonding is saturated, i.e. no more than some characteristic 
number of bonds N may be permitted per one molecule. The same features may be suggested 
for the bonds in polymers, so the statistical description of polymer solutions may be applied 
to the molecular liquids with bonding. The difference is that in the theory of polymers, one 
supposes the polymer chain very long, while in our case, the chain of bonded molecules may have 
an arbitrary length. 

The idea of the model arises from the works of Baxter [22], Nikomarov and Obukhov 
and may be formulated as follows. 

Consider spatial net of sites which approximately corresponds to the spacing of molecules in 
the liquid. Neighboring sites are connected with each other by edges. The number of nearest 
neighbors we denote as 7. Let us numerate the sites by Greece indexes, and the edges by Latin 
ones. Each site a is characterized by variable n a , which is the filling number 



Also, real scalar variable ipi corresponds to the edge numbered by %. Let us consider the partition 



2 The model. 




1, if the site a is occupied by the molecule 
0, if the site a is empty 





(2) 
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where the summation and integration goes over all possible configurations {n, ip} of variables 
n a ,ipi. The "effective Hamiltonian" F{n, ip} is defined by 



F{n, 0} = -- ri a J a pn p - | ^ n a + -J^ K$ - ^ In [1 + n a R a {^)\ (3) 

a,/3 a i a 

R a (ip) = ai ^ + ° 2 ^£ + • • • + aN Yl ^a^S ■ ■ ■ V'iJT- ( 4 ) 

t a L a/ t a t a7- t a 

Here, indexes numerate the edges which are connected with the site a. Spatial 
matrix J a p is 

T { J/T, if a,j3 are the nearest neighbors , , 

af3 ~ \ 0, otherwise ' [ } 

where J is the energy of attractive non - directed (Van der Waals) interaction between molecules, 
and fi is their chemical potential. Thus, the first and the second terms in (JHJ) correspond to the 
lattice - gas model [22] and allow one to describe the density of the system. The value K is 

K = exp(U/T), (6) 

where U is the energy of quasi - chemical bond. In (0 EI), T = A^T^, where Tk,ks are the 
temperature and the Boltzmann constant respectively. The last two terms inQ allow one to 
describe bonding. To demonstrate this fact, let us rewrite (0) 



{n} [ i 



I (1 + n a R a (i>)) 



a/3 a 

and consider the functional integral over {ijj} at some given configuration of filling numbers {n}. 
Releasing the brackets in the product, one gets the sum over all possible series of products: 

J (1 + n a R a (ip)) = 1 + niai^ + niai^ 2l + • • • n^ajip^^ + (8) 

a 

Then, the integration with the Gaussian weight exp [— |if -0?] , produces all possible couplings 
between pairs of neighboring sites. The coupling takes place when the mathematical power of x/i 
on corresponding edge is two. Zero power corresponds to the edge without coupling. The ratio 
of the weight of coupled edge to the weight of uncoupled one is K~ x = exp(— £//T). Thus, the 
partition (JHJ) generates the sum over all possible configurations of occupied sites, connected by all 
possible couplings between nearest neighbors, which correspond to intermolecular bonds. Each 
site with k bonds along its edges has an additional weight a^- No more than one bond per edge, 
and no more than N bonds per site are permitted due to the structure of polynomial R a (ip), see 
®. 

Besides, each configuration of filling numbers {n} has the standard weight of the lattice - gas 
theory, exp[| J2 a/3 n a Japrip — §J2 a n a ]. Thus, in the frames of formalism considered, one can 
easily model an arbitrary molecular liquid by an appropriate choice of model parameters. Those 
are: 
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N 
J 
/' 



The energy of quasi - chemical bond 

The maximal number of bonds per molecule 

The energy of non - directed (Van der Waals) interaction 

The chemical potential of molecules 

The weight of m - bonded molecule 



Instead of filling numbers, it is more useful to deal with scalar variables without any limi- 
tations imposed. To do so, Hubbard - Stratonovich transform to conjugated field may be used. 
Then, the summation over filling numbers leads to 



Z = J Dil>D<pexp[-F{(p,ij)}] 

i a/3 a 



(9) 



To analyze functional integral one can apply traditional method: to find most probable 
configuration which minimizes the F{(p,ip}, then to investigate fluctuations around it, etc. In 
the next section, we describe the first step, which is the mean - field approximation. 



3 Mean - field approximation. 



The mean - field configuration < tp a >, < ^ > provides the minima of (jOJ and obeys following 
equations: 

OF dF 

&fc !<*«>.<*>- °' !<*>«>.<*>- °- ( 10 ) 

Instead of <p a , it is better to use the variable 

q a = E J «/3V/3, 
P 

because its mean - field value, 

w a =<q a >, (11) 

coincides with the mean value of filling number < n a > and should be understood as the mean 
concentration of the molecules. 

For homogeneous system, < ^ >= w a = w, equations (fTUjl may be written in explicit 
form 



exp[^ W + ^ + ln(l + i^))] 
l + expfe; + £ + ln(l + i?(^))] 



w = —j 12 



n s u R'(V) 

* = —e-?w- K —^- 13 

n e 1 + R(i&) v ' 

R{V) = a l7 ^ + a 2 7(T ~ 1 ^ 2 + ... + a N C^ N (14) 

where n s /n e is the ratio of the number of sites to the number of edges in the lattice and 7 is 
the number of nearest neighbors. Solutions of system (jl2H14j) determine the equilibrium values 
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of ty,w. For nonequilibrium ^f,w, expression gives the density of thermodynamic potential 
in the Landau theory 



He 

2n s 



■erf 



2T 



w — In 



1 + e^ 1 



(15) 



Solutions of system (|12H14|) correspond to different phases of the system. The stable phases 
correspond to minimums of (II 5|) . 

The first of equations (fT2"J) is exactly the equation which arises in the lattice - gas theory of 
critical point: 

Jl ■ - - 11 ■ — ■ R (n- (i6) 



In 



-=riv + B, B 



|+ln(l 



1 - w T 

Depending on parameters T, 5, this equation may have from one to three solutions on the interval 
(0, 1). The solution with w ~ 1 corresponds to the condensed phase (liquid), while the solution 
w ~ represents the gas. We are interested in the condensed phase behavior, far from the liquid 
- gas critical point. Thus, we consider the case when the equation (fT6|) has only one solution, 
w ~ 1. This may always be provided by an appropriate choice of parameter fi in fTT)|) . If the 
equation (jlfi)) has only one solution, then one gets the Landau theory with single order parameter 
\l/ and with thermodynamic potential 



/(*) 



2b 



-w (fy) - In 



l + e «»(*)+$ (! + #(¥)) 



(17) 



where w(^) is given by (fT7)|) . 6 = ^e~^, a = 4?. Function f(^) may have several minima, 
which corresponds to different local arrangements of bonds. These minima are determined by 
equation (fT^j) : 

R '(^) 

* = teW rrw (18) 

The formalism considered allows to describe a wide variety of intermolecular bonding in a 
liquid. The bonding has been modelled by the choice of the polynomial under the logarithm in 



R(q) = « l7 ^ + a 2 — — * 2 



N,t,N 



(19) 



The number N is the maximal number of bonds per molecule, and the coefficient a m is the weight 
of m-bonded molecule. Note, that P(^) is positive for most physical cases. Since the coupling 
arises due to the charge redistribution in the molecule, the maximal number of bonds is even. 
The odd number of bonds means that the charge distribution in the molecule is not symmetric, 
so 

«2n > Ct2n-1- 

It can be easily shown, that this provides positive values of R{^) at any 
At low temperatures, w(^>) ~ 1, and (|18|) may be rewritten as 



B 



exp 



U 
T 



i + R(yy 

Its graphic solution is presented on fig. la. The right side of the equation, 
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with some intermediate behavior, as it presented on fig. la. Since U < 0, the left side is a line 
which slope lowers with temperature lowering. Thus, if the temperature is low enough, there 
always exist nontrivial set of solutions, which starts from the minimum of f(^f) at large negative 
and completes with another minimum at large positive 

Here, some comments should be done. The mean - field approximation does not answer the 
question what is the physical difference between these two low temperature structures. To do 
so, one needs higher approximations. Similar situation takes place for the Ising model in two 
dimensions [213 • ( m 2D, the Ising model may be formulated in terms of scalar field which corre- 
sponds to magnetization. At the same time, it may be understood as a statistics of self - closed 
contours which divide areas with inverted spins. Mean magnetization gives poor information 
about contours). Nevertheless, some assertions may be done. In our case, absolute value of ^ is 
proportional to the mean length of quasi - polymers [T^]. The structure \& < is metastable, i.e. 
/(\& < 0) > /(\& > 0) (fig. lb). It may be stabilized by external field s conjugated to ip, i.e. by 
adding the term sip to the Hamiltonian. Such a term produces spatial points, where the quasi - 
polymers have been terminated. Thus, the difference between \P > and ^ < structures lies 
in the statistics of polymer tails. For the Bethe lattice, the metastable state does not exist (see 
the next section). Since for the Bethe lattice closed paths are prohibited, then it is natural to 
suggest that the metastable state may be characterized by a set of self - closed quasi - polymers. 

When the maximal number of bonds per molecule N > 4, additional solutions of equation (|18|) 
arise. Mathematically, these solutions arise due to kinks in the right hand side of IflHJ) . provided 
by high powers in the polynomial R{^)- As an example, the graphic solution of eq. ()18j) at 
some certain model parameters is presented at fig.lc, together with corresponding F(ty) (fig. Id) 
dependence. As can be easily understood, the existence of such solutions may result in the liquid 
- liquid phase transition. The model predicts possible phase diagrams in the [i — T plane, which 
are presented at fig.2. Phase diagrams of that type were suggested earlier [22l|23j. Note, that 
the diagram at fig. 2b demonstrates an interesting possibility for the system behavior: there is no 
first order transition in density parameter w, and \1/ occurs to be the relevant order parameter 
for the liquid - gas critical point. 

Thus, the mean - field approximation predicts two non - trivial features. First, there exists 
metastable state with negative This state appears at rather low temperatures. Second, high 
powers of R{^) which correspond to the case N > 4, may provide liquid - liquid phase transition 
at some sets of weights a n . The transition is the first order phase transition, which line terminates 
in the critical point at the [i — T plane. At this line, order parameter jumps from ^ to In 
the critical point, the jump vanishes. Since the absolute value of \1/ is proportional to the mean 
length of quasi - polymers, then higher \l/ corresponds to the higher number of active bonds. The 
density difference is small, but higher \l/ corresponds to the higher density also. 

To verify these results, in the next section we analyzed the exact solution of the model on 
the Bethe lattice. 

4 The model on the Bethe lattice. 

In this section, we put uo a = 1 for simplicity. All the results obtained may be easily derived for 
arbitrary concentration, but we are interested mostly in the high density limit. Thus, we deal 
with the partition 




a 



(21) 
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Figure 1: The graphic solution of equation (fTHJ). i.e. left and right sides of this equation, and 
corresponding plot of thermodynamic potential (JT7j) at T = T/U = 0.5 (a - graphic solution, 
b - thermodynamic potential) and T = T/U = 1.019 (c - graphic solution, d - thermodynamic 
potential). The calculations are done in reduced variables: n e /n s — A = 3, 7 = 6, a\ — . . . — 
(I4 = a = 0.015, J = J/U = 0.2, jl = fi/U = 1. Parameter iV = 4 in both cases. 

on the Bethe lattice. Here, 

i % 

In i^O; we used the source term with Si = s. Partition (|21jl may be treated as the generating 
function for mean values 

(22) 

For each edge % let us define function Q(ipi), which is the partition (|21|) for sublattice, obtained 
from the initial Bethe lattice by cut of the edge %. In the thermodynamic limit, the partition for 
initial lattice depends no on the choice of the edge and is given by 

Z = j # 4 e"f tf+^Q 2 ^). (23) 

Let us introduce the values 

A = J d^e-^ 2+s ^Q(^, (24) 
7 
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Figure 2: The jx — T phase diagrams at J = 0.7 (a) and J = 0.3, (b). Other parameters are 
n e /n s = A = 3, 7 = 6, a± — . . . — 04 = a — 0.015, N = 4. Liquid - liquid and liquid - gas 
critical points are denoted as Ci and C2 respectively 



B 



(25) 



At s = 0, these values are proportional to the probabilities of the edge to be with or without 
bond, respectively. Obviously, 



Qty) = F 1 (A,Bty + F 2 {A,B), 



(26) 



where explicit form of functions Fi(A, B) and F 2 (A, B) depends on model details (the number of 
nearest neighbors and the set of coefficients a n ). Integrating (|2fij) . one gets system of equations 
on values A, B: 



A = F 1 (A, B)J 2 (K, s) + F 2 (A, B)J X {K, s) 
B = F X (A, B)J X {K, s) + F 2 (A, B)J (K, s). 



Here, 



Jn(K, S) 



dipe 



cP" 1 



ds 



n-l 



se 2K 



From (j22jl . (Pjl . (|2Bjl . (|27jl one obtains 

2F X F 2 



, 2X 3F? + F*K 



FlK + F^K 2 ' 



F 2 + F 2 2 K J 
At s = system (|27|) takes the form 

A = F 1 (A,B)J 2 (K, 0) = Fi(A, B)\ph\K~ 



B = F 2 (A, B)J (K, 0) = F 2 (A, B)y/2nK-?, 



(27) 



(28) 



(29) 
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and 

2AB 3KA 2 + B 2 . , 

W = ^ + ffl ' &) = K*A* + KB* - (30) 

As can be seen from (I29J) . any mean value may be presented as a function of the parameter 
x = A/B: 

w = *Jrr <^ = 1S^f < 31 > 

The equation on x is 

Kx = y(x), (32) 

where y(x) = Fi(x)/F 2 (x). 

Any thermodynamic function may be reconstructed using (ip 2 ). For example, 

Using known relation for internal energy, 

£ = r 2 ^ (33) 

one gets 

E=-^K{^ 2 ). (34) 

Inserting F = — T In Z into (|33|) and integrating, one gets the free energy thermodynamic poten- 
tial 

T 

F(T) - F(T ) = || rfr^ <^(r)>, (35) 

where (^ 2 (r)) is defined by (gQ and (g2D - 

Thus, in the thermodynamic limit, any thermodynamic characteristic may be calculated using 
y(x). Now, let us determine this function for arbitrary model parameters. Those are the maximal 
number of bonds per molecule N, the number of nearest neighbors 7, and the set of coefficients 
a n . The only limitation imposed is 7 > N. Consider some site a and suggest that the integration 
in (JHJ) is done over all ip except those which correspond to the site a. If one of these edges is cut 
(for definition, the edge ipi), then initial lattice occurs to be divided in two sublattices (fig. 3). 

Then, Q(ipi) may be written as 



j (jJd^ie-^Q^i)] [1 + ai(V>i + • • • + V 7 ) + ^2(^1^2 + • • • + ^7-1^7) + ■ ■ ■] 



(36) 

In the thermodynamic limit, Q{ipi) depends no on i. Integration together with (J24j) . IJ25|) and 
ff2"6l gives 

F ± (A, B) = ai B^~ l + a 2 ( 7 - l)AB^~ 2 + ... + ajvC^A^ 1 W~ N 
F 2 (A, B) = B 1 " 1 + 01(7 - l)AB^ 2 + ... + a N Cl~ l A N B^ N ~\ (37) 
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Figure 3: The cut of the initial Bethe lattice across the edge leads to the equation (|3l)|) 



where C 7 ^ = n\/m\(n — m)\ - combinatoric coefficients. Since x = A/B, y(x) = Fi(x)/F 2 (x), 
then 

N-l 

y(x) = ^ — • (38) 

Q"niCm X m 

m=0 

Thus, equation (|3*2"J) is algebraic with power N—l. Qualitatively, its structure coincides with mean 
- field equation ()18j) . It should be mentioned, however, that (}3*2j) may be treated as a recurrence 
relation Kx n = y(x n _i), which connects x values for n and n — l levels of the Cayley tree, taken 
in the thermodynamic limit. Since any initial x = A/B is positive, then all stationary points of 
this recurrence are positive too. Thus, on the Bethe lattice, the low - temperature metastable 
state does not exist. 



5 Conclusion 

Here, we point out some general results and features which make the model rather attractive. 

1. In mean field approximation, the metastable low - temperature state is predicted. This 
state is universal, i.e. appears for any set of model parameters. It becomes stable only 
at zero temperature, but may be stabilized by an appropriate external field. For example, 
phase separation in two phases with positive and negative \l/ may be achieved by impurity 
addition. Recently, similar separation was observed experimentally [21]. For the Bethe 
lattice, the metastable state does not exist. Since for the Bethe lattice closed paths are 
prohibited, then it is natural to suggest that the metastable state may be characterized by 
a set of self - closed quasi - polymers. 

2. Liquid - liquid phase transition between phases with different density of bonds is predicted 
both in mean - field approximation and on the Bethe lattice. Possible phase diagrams 
correspond to those predicted by computer simulations [221 EHJ The phase transition is 



10 



the first order phase transition which takes place at some special sets of model parameters. 
Also, the set of model parameters may provide not phase transition but sharp crossover 
from one mean density of bonds to another. In that case, the order parameter temperature 
dependence demonstrates a kink. Thus, all \I/ - dependent properties should demonstrate 
a kink too. Such a kinks have often been observed in molecular liquids [7]-|14j. 

3. The model considered allows one to describe an arbitrary set of molecules and their bonding. 
For example, different conformations of bonds may be easily accounted by an appropriate 
choice of R(ip). 

4. Actually, the model considered is the vertex model. Such models have often been used in 
statistical physics. For example, Ising model may be formulated as the vertex model |25j . 
Thus, the methods to investigate the model behavior are developed rather well. 

The work is supported by RFBR (grants N 04-02-96095, 06-02-17269). 
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